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Search for the associated production of the Higgs boson with a top quark pair in 
multilepton final states with the ATLAS detector 


The ATLAS Collaboration 


Abstract 

A search for the associated production of the Higgs boson with a top quark pair is performed in multilepton final states using 
20.3 fb 1 of proton-proton collision data recorded by the ATLAS experiment at yfs - 8 TeV at the Large Hadron Collider. Five 
final states, targeting the decays H —> WW*, tt, and ZZ*, are examined for the presence of the Standard Model (SM) Higgs 
boson: two same-charge light leptons (e or /./) without a hadronically decaying r lepton; three light leptons; two same-charge light 
leptons with a hadronically decaying r lepton; four light leptons; and one light lepton and two hadronically decaying r leptons. 
No significant excess of events is observed above the background expectation. The best fit for the ttH production cross section, 
assuming a Higgs boson mass of 125 GeV, is 2.1+}^ times the SM expectation, and the observed (expected) upper limit at the 95% 
confidence level is 4.7 (2.4) times the SM rate. The p-value for compatibility with the background-only hypothesis is 1.8cr; the 
expectation in the presence of a Standard Model signal is 0.9cr. 


1. Introduction 

The discovery of a new particle H with a mass of about 
125 GeV in searches for the Standard Model (SM) [1-3] Higgs 
boson [4-7] at the LHC was reported by the ATLAS [8] and 
CMS [9] Collaborations in July 2012. The particle has been 
observed in the decays H —> yy [ 10, 11], H —> ZZ* —> 4f 
[12, 13], and H — > WW* —* Mv [14, 15], and evidence has 
been reported for H — » tt [16, 17], consistent with the rates 
expected for the SM Higgs boson. 

The observation of the process in which the Higgs boson is 
produced in association with a pair of top quarks (ttH) would 
permit a direct measurement of the top quark-Higgs boson 
Yukawa coupling in a process that is tree-level at the lowest 
order, which is otherwise accessible primarily through loop ef¬ 
fects. Having both the tree- and loop-level measurements would 
allow disambiguation of new physics effects that could affect 
the two differently, such as dimension-six operators contribut¬ 
ing to the ggH vertex. This letter describes a search for the SM 
Higgs boson in the ttH production mode in multilepton final 
states. The five final states considered are: two same-charge- 
sign light leptons (e or p) with no additional hadronically de¬ 
caying r lepton; three light leptons; two same-sign light leptons 
with one hadronically decaying r lepton; four light leptons; and 
one light lepton with two hadronically decaying r candidates. 
These channels are sensitive to the Higgs decays H —» WW*, 
rr, and ZZ* produced in association with a top quark pair de¬ 
caying to one or two leptons. A similar search has been per¬ 
formed by the CMS Collaboration [18]. 

The selections of this search are designed to avoid overlap 
with ATLAS searches for ttH in H —> bb [19] and H —> yy [20] 
decays. The main backgrounds to the signal arise from ft pro¬ 
duction with additional jets and non-prompt leptons, associated 
production of a top quark pair and a vector boson W or Z (col¬ 


lectively denoted ttV), and other processes where the electron 
charge is incorrectly measured or where quark or gluon jets are 
incorrectly identified as t candidates. 

2. ATLAS detector and dataset 

The features of the ATLAS detector [21] most relevant to 
this analysis are briefly summarized here. The detector consists 
of an inner tracking detector system surrounded by a supercon¬ 
ducting solenoid, electromagnetic and hadronic calorimeters, 
and a muon spectrometer. Charged particles in the pseudorapid¬ 
ity 1 range \r/\ < 2.5 are reconstructed with the inner tracking 
detector, which is immersed in a 2 T magnetic field parallel to 
the detector axis and consists of pixel and strip semiconductor 
detectors as well as a straw-tube transition radiation tracker. 
The solenoid is surrounded by a calorimeter system covering 
\rj\ < 4.9, which provides three-dimensional reconstruction of 
particle showers. Lead/liquid-argon (LAr) sampling technology 
is used for the electromagnetic component. Iron/scintillator- 
tile sampling calorimeters are used for the hadronic component 
for \r/\ < 1.7, and copper/LAr and tungsten/LAr technology is 
used for I 77 I > 1.5. Outside the calorimeter system, air-core 
toroids provide a magnetic field for the muon spectrometer. 
Three stations of precision drift tubes and cathode-strip cham¬ 
bers provide a measurement of the muon track position and 

1 The ATLAS experiment uses a right-handed coordinate system with its 
origin at the nominal interaction point (IP) in the centre of the detector, and 
the z-axis along the beam line. The x-axis points from the IP to the centre 
of the LHC ring, and the y-axis points upwards. Cylindrical coordinates (r, 0) 
are used in the transverse plane, 0 being the azimuthal angle around the z- 
axis. Observables labelled “transverse” are projected onto the x—y plane. The 
pseudorapidity is defined in terms of the polar angle 6 as 77 = - In tan 6/2. The 
transverse momentum is defined as pj = p sin 6 = pi cosh 77 , and the transverse 
energy Ej has an analogous definition. 
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curvature in the region \r/\ < 2.1. Resistive-plate and thin-gap 
chambers provide muon triggering capability up to \q\ = 2.4. 

This search uses data collected by the ATLAS experiment in 
2012 at a centre-of-mass energy of Vs = 8 TeV. All events con¬ 
sidered were recorded while the detector and trigger systems 
were fully functional; the integrated luminosity of this dataset 
is 20.3 fb -1 . 

3. Cross sections for signal and background processes 

The cross section for the production of ttH in pp collisions 
has been calculated at next-to-leading order (NLO) in quantum 
chromodynamics (QCD) [22-26]. Uncertainties on the cross 
section are evaluated by varying the renormalization and fac¬ 
torization scales by factors of two and by varying the input par- 
ton distribution functions (PDF) of the proton. A Higgs boson 
mass of inn = 125 GeV is assumed; this gives a predicted ttH 
production cross section at — 8 TeV of 129 (scale) +10 
(PDF) fb [27], This assumed Higgs boson mass is consistent 
with the combined ATLAS and CMS measurement [28], 

In this letter the associated production of single top quarks 
with a Higgs boson is considered a background process and set 
to the Standard Model rate. The production of tHqb and t/l W 
is taken into account. In the Standard Model these rates are 
very small compared to ttH production. These processes are 
simulated with the same parameters as used by the ATLAS ttH, 
H —> yy search [20], The cross sections for both are computed 
using the MG5 _aMC@NLO generator [29] at NLO in QCD. 
For tHqb, the renormalization and factorization scales are set to 
75 GeV and the process is computed in the four-flavour scheme, 
yielding cr(jHqb) =11.2 ^ (scale) (PDF) fb. For tHW, 
dynamic factorization and renormalization scales are used, and 
the process is computed in the five-flavour scheme; the result 
is cr(tHW) = 4.7 (scale) (PDF) fb. The interference 
of tHW production with ttH, which appears at NLO for tHW 
in diagrams with an additional fc-quark in the final state, is not 
considered. 

The production of ttW and ttiZ/y 1 ') —> ttH H yield mul¬ 
tilepton final states with Z?-quarks and are major backgrounds to 
the ttH signal. For simplicity of notation the latter process is re¬ 
ferred to as ttZ throughout this letter with off-shell Z and photon 
components also included except where noted otherwise. The 
ttW process includes both ttW 1 and ttW components. Next-to- 
leading-order cross sections are used for ttW [30] and ttZ [31]. 
The MG5 _aMC@NLO generator is used to reproduce the QCD 
scale uncertainties of these calculations and determine uncer¬ 
tainties due to the PDF. For ttW production the value 232 ± 28 
(scale) +18 (PDF) fb is used, and for ttZ production 2 the value 
is 206 ± 23 (scale) ± 18 (PDF) fb. 

The associated production of a single top quark and a Z bo¬ 
son is a subleading background for the most sensitive channels. 


2 The NLO cross section is only evaluated for ttZ production with on-shell 
Z. The cross section obtained for rf(Z/y*) production including off-shell Z/y* 
contributions in a leading-order simulation is scaled by a A”-factor of 1.35 ob¬ 
tained as the ratio of NLO and LO on-shell cross sections. The A'-factor differs 
from that of Ref. [31] due to a different choice of PDF. 


The cross section has been calculated at NLO for the t- and s- 
channels [32], The resulting values used in this work are 160 
± 7 (scale) ± 11 (PDF) fb for tZ and 76 + 4 (scale) ± 5 (PDF) 
fb for tZ. The cross section for the production of tWZ is com¬ 
puted at leading order (LO) using the MadGraph v5 generator 
[33] and found to be 4.1 fb. 

The cross section for inclusive production of vector boson 
pairs WW, WZ, and ZZ is computed using MCFM [34], Con¬ 
tributions from virtual photons and off-shell Z bosons are in¬ 
cluded. The uncertainties on the acceptance for these processes 
in the signal regions (which favour production with additional 
b- or c-quarks) dominate over the inclusive cross-section un¬ 
certainty (see Section 7.2) and so the latter is neglected in the 
analysis. 

The inclusive tt cross section is calculated at next-to-next- 
to-leading order (NNLO) in QCD which includes resummation 
of next-to-next-to-leading logarithmic (NNLL) soft gluon terms 
using Top++ [35], yielding 253+[j pb for yfs = 8 TeV. The 
single-top-quark samples are normalized to the approximate 
NNLO theoretical cross sections [36-38] using the MSTW2008 
[39] NNLO PDF set. The production of Z —> CC+jets and 
W —> (V+jets is normalized using NNLO cross sections as com¬ 
puted by FEWZ [40], 

4. Event generation 

The event generator configurations used for simulating the 
signal and main background processes are shown in Table 1. 
Additional information is given below. 

The ttH signal event simulation samples contain all Higgs 
boson decays with branching fractions set to values computed 
at NNLO in QCD [26, 66-69], The factorization (pp) and renor¬ 
malization (pH scales are set to m, + m H jl. Higgs boson and 
top quark masses of 125 and 172.5 GeV, respectively, are used. 
These samples are the same as those used by other ATLAS ttH 
searches [19, 20]. 

Production of single top quarks with Higgs bosons is sim¬ 
ulated as follows. For tHqb, events are generated at leading 
order with MadGraph in the four-flavour scheme. For tHW, 
events are generated at NLO with MG5_aMC@NLO in the 
five-flavour scheme. Higgs boson and top quark masses are set 
as for ttH production. 

The main irreducible backgrounds are production of ttW 
and ttZ (ttV). For the ttW process, events are generated at lead¬ 
ing order with zero, one, or two extra partons in the final state, 
while for ttZ zero or one extra parton is generated. The import¬ 
ant contribution from off-shell y*/Z —> HH is included. The tZ 
process is simulated with the same setup, without extra partons. 

For diboson processes, the full matrix element for Ht pro¬ 
duction, including y* and off-shell Z contributions, is used. The 
Sherpa qq and qg samples include diagrams with additional par- 
tons in the final state at the matrix-element (ME) level, and in¬ 
clude b- and c-quark mass effects. Sherpa was found to have 
better agreement with data than Powheg for WZ, while the 
Sherpa and Powheg descriptions of ZZ production are similar. 

A ff+jets sample generated with the Powheg NLO gener¬ 
ator [61] is used; the top quark mass is set to 172.5 GeV. Small 
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Table 1: Configurations used for event generation of signal and background processes. If only one parton distribution function is 
shown, the same one is used for both the matrix element (ME) and parton shower generators; if two are shown, the first is used for 
the matrix element calculation and the second for the parton shower. “Tune” refers to the underlying-event tune of the parton shower 
generator. “Pythia 6 ” refers to version 6.425; “Pythia 8 ” refers to version 8.1; “Herwig++” refers to version 2.6; “MadGraph” 
refers to version 5; “Alpgen” refers to version 2.14; “Sherpa” refers to version 1.4; “gg2ZZ” refers to version 2.0. 


Process 

ME Generator 

Parton Shower 

PDF 

Tune 

ttH 

HELAC-Oneloop [41, 42] 

+ Powheg-BOX [48-50] 

Pythia 8 [43] 

CT10 [44J/CTEQ6L1 [45, 46] 

AU2 [47] 

tHqb 

MadGraph [33] 

Pythia 8 

CT10 

AU2 

tHW 

MG5_aMC@NLO [29] 

Herwig++ [51] 

CT10/MRST LO** [52] 

UE-EE-4 [53] 

ttW + <2 partons 

MadGraph 

Pythia 6 [54] 

CTEQ6L1 

AUET2B [55] 

ff(Z/y*) + < 1 parton 

MadGraph 

Pythia 6 

CTEQ6L1 

AUET2B 

Kz/r ) 

MadGraph 

Pythia 6 

CTEQ6L1 

AUET2B 

qq, qg -> WW, WZ 

Sherpa [56] 

Sherpa 

CT10 

Sherpa default 

qq —> qqWW, qqWZ, qqZZ 

Sherpa 

Sherpa 

CT10 

Sherpa default 

qq, qg -> ZZ 

Powheg-BOX [57] 

Pythia 8 

CT10 

AU2 

gg -> zz 

gg2ZZ [58] 

Herwig [59] 

CT10 

AUET2 [60] 

tt 

Powheg-BOX [61] 

Pythia 6 

CT10/CTEQ6L1 

Perugia2011C [62] 

s-, t-channel, Wt single top 

Powheg-BOX [63, 64] 

Pythia 6 

CT10/CTEQ6L1 

Perugia2011C 

Z —» TT + < 5 partons 

Alpgen [65] 

Pythia 6 

CTEQ6L1 

Perugia2011C 

W —> £v+ < 5 partons 

Alpgen 

Pythia 6 

CTEQ6L1 

Perugia2011C 


corrections to the tt system and top quark p\ spectra are ap¬ 
plied based on discrepancies in differential distributions ob¬ 
served between data and simulation at 7 TeV [70], Double¬ 
counting between the tt and VV7 single top production final states 
is eliminated using the diagram-removal method [71], 

Samples of Z —> TT^+jets and W —> (V+jets events are 
generated with up to five additional partons using the Alp¬ 
gen v2.14 [65] leading order (LO) generator. Samples are 
merged with matrix element-parton shower overlaps removed 
using MLM matching [72]. Production of b- and c-quarks is 
also computed at matrix-element level, and overlaps between 
ME and parton shower production are handled by separating 
the kinematic regimes based on the angular separation of addi¬ 
tional heavy partons. The resulting “light” and “heavy” flavour 
samples are normalized by comparing the resulting ^-tagged jet 
spectra with data. 

All simulated samples with Pythia 6 and Herwig [59] par- 
ton showering use Photos 2.15 [73] to model photon radiation 
and Tauola 1.20 [74] for r decays. The Herwig++ samples 
model photon radiation with Photos but use the internal r decay 
model. Samples using Pythia 8.1 and Sherpa use those gener¬ 
ators’ internal r lepton decay and photon radiation generators. 
For Herwig samples, multiple parton interactions are modelled 
with Jimmy [75], 

Showered and hadronized events are passed through sim¬ 
ulations of the ATLAS detector (either full GEANT4 [76] 
simulation or a hybrid simulation with parameterized calori¬ 
meter showers and GEANT4 simulation of the tracking sys¬ 
tems [77, 78]). Additional minimum-bias pp interactions 
(pileup) are modelled with the Pythia 8.1 generator with the 
MSTW2008 LO PDF set and the A2 tune [79], They are added 
to the signal and background simulated events according to the 
luminosity profile of the recorded data, with additional overall 


scaling to achieve a good match to observed calorimetry and 
tracking variables. The contributions from pileup interactions 
both within the same bunch crossing as the hard-scattering pro¬ 
cess and in neighbouring bunch crossings are included in the 
simulation. 

5. Object selection 

Electron candidates are reconstructed from energy clusters 
in the electromagnetic calorimeter associated with reconstruc¬ 
ted tracks in the inner detector. They are required to have 
l^ciusterl < 2.47. Candidates in the transition region 1.37 < 

I ^cluster I < 1-52 between sections of the electromagnetic calor¬ 
imeter are excluded. A multivariate discriminant based on 
shower shape and track information is used to distinguish elec¬ 
trons from hadronic showers [80, 81]. Only electron candid¬ 
ates with transverse energy E T greater than 10 GeV are con¬ 
sidered. To reduce the background from non-prompt electrons, 
i.e. from decays of hadrons (including heavy flavour) produced 
in jets, electron candidates are required to be isolated. Two 
isolation variables, based on calorimetric and tracking vari¬ 
ables, are computed. The first (£’“ ne ) is based on the sum of 
transverse energies of calorimeter cells within a cone of radius 
A R = -y/(A tf>) 2 + (A/ 7) 2 = 0.2 around the electron candidate 
direction. This energy sum excludes cells associated with the 
electron and is corrected for leakage from the electromagnetic 
shower and ambient energy in the event. The second (/tj onc ) is 
defined based on tracks with pj > 1 GeV within a cone of radius 
A R = 0.2 around the electron candidate. Both isolation energies 
are separately required to be less than 0.05 x Ej. The longitud¬ 
inal impact parameter of the electron track with respect to the 
selected event primary vertex, multiplied by the sine of the polar 
angle, |zo sin G\, is required to be less than 1 mm. The transverse 
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impact parameter divided by the estimated uncertainty on its 
measurement, \da\lcr{da), must be less than 4. If two electrons 
closer than AR = 0.1 are selected, only the one with the higher 
p-y is considered. An electron is rejected if, after passing all the 
above selections, it lies within A R - 0.1 of a selected muon. 

Muon candidates are reconstructed by combining inner de¬ 
tector tracks with track segments or full tracks in the muon 
spectrometer [82]. Only candidates with \q\ < 2.5 and p T > 
10 GeV are kept. Additionally, muons are required to be sep¬ 
arated by at least AR > 0.04 + (10 GeV)//? T>; , from any selec¬ 
ted jets (see below for details on jet reconstruction and selec¬ 
tion). The cut value is optimized to maximize the acceptance 
for prompt muons at a fixed rejection factor for non-prompt and 
fake muon candidates. Furthermore, muons must satisfy sim¬ 
ilar Ej° ne and /?™ ne isolation criteria as for electrons, with both 
required to be less than 0.10 x pj. The value of |zo sin 0\ is re¬ 
quired to be less than 1 mm, while \do\/cr(do) must be less than 
3. 

Hadronically decaying t candidates (rh a d) are reconstruc¬ 
ted using clusters in the electromagnetic and hadronic calori¬ 
meters. The r candidates are required to have pj greater than 
25 GeV and //| < 2.47. The number of charged tracks as¬ 
sociated with the r candidates is required to be one or three 
and the charge of the r candidates, determined from the associ¬ 
ated tracks, must be +1. The r identification uses calorimeter 
cluster and tracking-based variables, combined using a boos¬ 
ted decision tree (BDT) [83]. An additional BDT which uses 
combined calorimeter and track quantities is employed to reject 
electrons reconstructed as one-prong hadronically decaying r 
leptons. 

Jets are reconstructed from calibrated topological 
clusters [21] built from energy deposits in the calorimet¬ 
ers, using the anti-k, algorithm [84-86] with a radius parameter 
R = 0.4. Prior to jet finding, a local cluster calibration 
scheme [87, 88] is applied to correct the topological cluster 
energies for the effects of non-compensating calorimeter 
response, inactive material and out-of-cluster leakage. The 
jets are calibrated using energy and //-dependent calibration 
factors, derived from simulations, to the mean energy of stable 
particles inside the jets. Additional corrections to account for 
the difference between simulation and data are derived from 
in-situ techniques [89, 90], After energy calibration, jets are 
required to have pj > 25 GeV and \r/\ < 2.5. 

To reduce the contamination from jets originating in pp in¬ 
teractions within the same bunch crossing (pileup), the scalar 
sum of the p-y of tracks matched to the jet and originating from 
the primary vertex must be at least 50% of the scalar sum of 
the p-y of all tracks matched to the jet. This criterion is only 
applied to jets with p-y < 50 GeV (those most likely to originate 
from pileup) and \q\ < 2.4 (to avoid inefficiency at the edge of 
tracking acceptance). 

The calorimeter energy deposits from electrons are typically 
also reconstructed as jets; in order to eliminate double count¬ 
ing, any jets within AR = 0.3 of a selected electron are not 
considered. 

Jets containing //-hadrons are identified (//-tagged) via a 
multivariate discriminant [91] that combines information from 


the impact parameters of displaced tracks with topological 
properties of secondary and tertiary decay vertices reconstruc¬ 
ted within the jet. The working point used for this search cor¬ 
responds to approximately 70% efficiency to tag a //-hadron jet, 
with a light-jet mistag rate of ~ 1 % and a charm-jet rejection 
factor of 5, as determined for //-tagged jets with p-y of 20-100 
GeV and \rj\ < 2.5 in simulated 11 events. To avoid inefficien¬ 
cies associated with the edge of the tracking coverage, only jets 
with |t 7 [ < 2.4 are considered as possible //-tagged jets in this 
analysis. The efficiency and mistag rates of the //-tagging al¬ 
gorithm are measured in data [91, 92] and correction factors 
are applied to the simulated events. 

6. Event selection and classification 

All events considered in this analysis are required to pass 
single-lepton (e or p) triggers. These achieve their maximal 
plateau efficiency for lepton p-y > 25 GeV. 

This analysis primarily targets the H —> WW* and tt decay 
modes. Considering the decay of the tt system as well, these 
ttH events contain either WWWWbb or rrWWbb. The strategy 
is to target final states that cannot be produced in tt decay alone 
— i.e., three or more leptons, or two same-sign leptons — thus 
suppressing what would otherwise be the largest single back¬ 
ground. 

The analysis categories are classified by the number of 
light leptons and hadronic r decay candidates. The leptons 
are selected using the criteria described earlier. Events are ini¬ 
tially classified by counting the number of light leptons with 
pj > 10 GeV. At least one light lepton is required to match a 
lepton selected by the trigger system. After initial sorting into 
analysis categories, in some cases the lepton selection criteria 
are tightened by raising the pj threshold, tightening isolation 
selections or restricting the allowed I/ 7 I range, as explained in 
the following per-category descriptions. The analysis includes 
five distinct categories: two same-sign light leptons with no Th a d 
(2/’0rh a d), three light leptons (3T), two same-sign light leptons 
and one Th a d (2A1 T| la d), four light leptons (4(j, and one light 
lepton and two Th a d (l/^Thad). The categories with Th a d candid¬ 
ates target the H —> tt decay; the others are primarily sensitive 
to H —> WW* with a very small contribution from H —> ZZ*. 
The contributions to each category from different Higgs boson 
decay modes are shown in Table 2. These selection criteria 
ensure that an event can only contribute to a single category. 
The contamination from gluon fusion, vector boson fusion, and 
associated VH production mechanisms for the Higgs boson is 
predicted to be negligible. Summed over all categories, the 
total expected number of reconstructed signal events assum¬ 
ing Standard Model ttH production is 10.2, corresponding to 
0.40% of all produced ttH events. The detailed criteria for each 
category are described below. 

6.1. 2 / 0 Th a d categories 

Selected events are required to include exactly two light 
leptons, which must have the same charge. Events with Th a d 
candidates are vetoed. To reduce the background from non¬ 
prompt leptons, the leading (subleading) lepton is required to 
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Table 2: Fraction of the expected ttH signal arising from dif¬ 
ferent Higgs boson decay modes in each analysis category. The 
six 2f0rhad categories are combined together, as are the two 41 
categories. The decays contributing to the “other” column are 
dominantly H —»pp and H —> bb. Rows may not add to 100% 
due to rounding. 



Higgs boson decay mode 

Category 

WW* 

TT 

ZZ* 

Other 

2F0r ha d 

80% 

15% 

3% 

2% 

3 1 

74% 

15% 

7% 

4% 

2^1 Thad 

35% 

62% 

2% 

1 % 

4 e 

69% 

14% 

14% 

4% 

lF2T ha d 

4% 

93% 

0% 

3% 


satisfy p T > 25 (20) GeV, and the muon isolation require¬ 
ments are tightened to E™ ne /p T < 0.05 and p™ ne /pj < 0.05. 
The angular acceptance of electron candidates is restricted to 
\rj\ < 1.37 in order to suppress tt background events where the 
sign of the electron charge is misreconstructed, as the charge 
misidentification rate increases at high pseudorapidity. 

In order to suppress the lower-multiplicity fF+jets and ttW 
backgrounds, events must include at least four reconstructed 
jets. In order to suppress diboson and single-boson back¬ 
grounds, at least one of these jets must be /^-tagged. The se¬ 
lected events are separated by lepton flavour (e ± e ± , e ± p ± , and 
p ± p ± ) and number of jets (exactly four jets, at least five jets) 
into six categories with different signal-to-background ratio, 
resulting in higher overall sensitivity to the ttH signal. 

6.2. 3£ category 

Selected events are required to include exactly three light 
leptons with total charge equal to ±1. Candidate events 
arising from non-prompt leptons overwhelmingly originate as 
opposite-sign dilepton events with one additional non-prompt 
lepton. As a result, the non-prompt lepton is generally one of 
the two leptons with the same charge. To reduce these back¬ 
grounds, a higher momentum threshold pj > 20 GeV is applied 
to the two leptons with the same charge. No requirements are 
imposed on the number of T haiJ candidates. In order to suppress 
the ff+jets and ttV backgrounds, selected events are required 
to include either at least four jets of which at least one must 
be /^-tagged, or exactly three jets of which at least two are b- 
tagged. To suppress the ttZ background, events that contain an 
opposite-sign same-flavour lepton pair with the dilepton invari¬ 
ant mass within 10 GeV of the Z mass are vetoed. Events con¬ 
taining an opposite-sign lepton pair with invariant mass below 
12 GeV are also removed to suppress background from reson¬ 
ances that decay to light leptons. 

6.3. 2Flrhad category 

Selected events are required to include exactly two light 
leptons, with the same charge and leading (subleading) pj > 
25 (15) GeV, and exactly one hadronic r candidate. The re¬ 
constructed charge of the Thad candidate has to be opposite 


to that of the light leptons. In order to reduce fF+jets and 
ttV backgrounds, events must include at least four reconstruc¬ 
ted jets. In order to suppress diboson and single-boson back¬ 
grounds, at least one jet must be /?-tagged. To suppress the 
Z —> F + F“+jets background, events with dielectron invariant 
mass within 10 GeV of the Z mass are vetoed. 

6.4. 4£ categories 

Selected events are required to include exactly four light 
leptons with total charge equal to zero and leading (subleading) 
pi > 25 (15) GeV. No requirements are applied on the number 
of Thad candidates. In order to suppress the fF+jets and ttV back¬ 
grounds, the selected events are required to include at least two 
jets of which at least one must be /?-tagged. To suppress the ttZ 
background, events that contain an opposite-sign same-flavour 
lepton pair with dilepton invariant mass within 10 GeV of the Z 
mass are vetoed. In order to suppress background contributions 
from resonances that decay to light leptons, all opposite-sign 
same-flavour lepton pairs are required to have a dilepton invari¬ 
ant mass greater than 10 GeV. The four-lepton invariant mass 
is required to be between 100 and 500 GeV, which gives high 
acceptance for ttH. H —> WW* —> tv tv. but rejects Z —> 4C 
and high-mass ttZ events. Selected events are separated by the 
presence or absence of a same-flavour, opposite-sign lepton pair 
into two categories, referred to respectively as the Z-enriched 
and Z-depleted categories. In both cases the Z mass veto is 
applied, but background events in the Z-enriched category can 
arise from off-shell Z and y* —* t + t~ processes while in the 
Z-depleted category these backgrounds are absent. 

6.5. lF2rhad category 

Selected events are required to include exactly one light 
lepton with pj > 25 GeV and exactly two hadronic r candid¬ 
ates. The Th a d candidates must have opposite charge. In order to 
suppress the fF+jets and ttV backgrounds, events must include 
at least three reconstructed jets. In order to suppress diboson 
and single-boson backgrounds, at least one of the jets must be 
/j-tagged. This final state is primarily sensitive to H —> r + r 
decays, allowing use of the invariant mass of the visible decay 
products of the Thadfhad system (m v ; s ) as a signal discriminant. 
Signal events are required to satisfy 60 < /;i vls < 120 GeV. 

7. Background estimation 

Important irreducible backgrounds include ttV and diboson 
production and are estimated from MC simulation. Validation 
regions enriched in these backgrounds are used to verify proper 
modelling of data by simulation. Reducible backgrounds are 
due to non-prompt lepton production and electron charge mis¬ 
identification, and are estimated from data, with input from sim¬ 
ulation in some categories. In the 1 F2r ha d category the primary 
concern is fake Th a d candidates, which are modelled using sim¬ 
ulation and validated against a data-driven estimate. 
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Figure 1: The spectrum of the number of jets expected and observed in the ttZ (left) and ttW (right) validation regions (VR). The 
hatched band represents the total uncertainty on the background prediction in each bin. The “non-prompt” backgrounds are those 
with a lepton arising from a hadron decay or from a photon conversion in detector material. Rare processes include tZ, ttWW, 
triboson, tttt, and tH production. The overlaid red line corresponds to the ttH signal predicted by the SM. 


7.1. ttV and tZ 

The primary backgrounds with prompt leptons stem from 
the production of ttW and ttZ. The ttW background tends to 
have lower jet multiplicity than the signal and so the leading 
contribution comes from events with additional high-px jets; 
it is the major ttV contribution in the 2fOT| lad categories and 
comparable to ttZ in the 2C\ T ha( | category. The ttZ process has 
similar multiplicity to the ttH signal but can only contribute 
to the signal categories when the Z boson decays leptonically, 
so the on-shell contribution can be removed by vetoing events 
with opposite-sign dilepton pairs with invariant mass near the 
Z pole. This is the larger of the two ttV contributions for the 
3C, AC, and 1 f'2T ivMi categories. The tZ process makes a sublead¬ 
ing contribution to both channels. A validation region is used 
to verify the modelling of ttZ using on-shell Z decays. Agree¬ 
ment is seen within the large statistical uncertainty. No region 
of equivalent purity and statistical power exists for ttW pro¬ 
duction; nevertheless the expectations are cross-checked with a 
validation region defined with the 2f0r ha d selection except with 
two or more Z?-tagged jets and either two or three jets, where the 
ttW purity is * 30%, and are found to be consistent within un¬ 
certainties. The spectra of the number of jets in these validation 
regions is shown in Fig. 1. 

Uncertainties on the ttV background contributions arise 
from both the overall cross section uncertainties (see Section 3) 
and the acceptance uncertainties. The latter are estimated by 
comparing particle-level samples after showering produced by 


three different pairs of generators: a) the nominal MadGraph 
LO merged sample versus an equivalent LO merged sample 
generated with Sherpa 2.1.1, to account for ME-parton shower 
matching effects; b) the LO merged Sherpa sample versus a 
Sherpa+OpenLoops [93] NLO sample, to compare LO merged 
and NLO acceptance; and c) MG5_aMC@NLO with Pythia 
8 parton shower versus Herwig++ parton shower, to compare 
px-ordered versus angular-ordered parton showers. Each of 
these variations is input independently into the final fit. When 
summed in quadrature they have an impact of 5-23% depend¬ 
ing on the category and background source (ttW versus ttZ). 
Uncertainties arising from changes in the acceptance due to the 
choice of QCD scale and PDL are also evaluated; these have an 
impact of 1.3-6.7% for scale and 0.9—4.8% for PDL 

7.2. Other prompt lepton contributions 

Other backgrounds with prompt leptons arise from multibo¬ 
son processes ( WZ , ZZ, and triboson production) in association 
with heavy-flavour jets, or with a misidentified light-flavour jet. 
The main process affecting the final result is WZ + jets. Val¬ 
idation regions with three leptons including a Z candidate and 
either zero or one Z?-tagged jet are studied. The number of jets 
in WZ + Ob events is reproduced well in the highly populated 
bins (up to 4 jets), leading to the conclusion that the jet ra¬ 
diation spectrum is well modelled. The dominant uncertainty 
on the prediction in the signal region is expected to arise from 
the WZ + b cross section. Data constrain this component with 
roughly 100% uncertainty. As a result a 100% uncertainty is 
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assigned to the WZ + b cross section, giving a 50% uncertainty 
on the total WZ yield, correlated across categories. The cross 
sections for production of WW + b and ZZ + b are also assigned 
50% uncertainties; these have negligible impact on the final res¬ 
ult. 

7.3. Charge sign misidentification 

The process <?* —» e ± y —* e ± e + e~ occurring in detector ma¬ 
terial can result in an electron produced with nearly the same 
momentum as the parent electron but with opposite charge. In 
these cases the observed electron has opposite charge to that 
of the primary electron (charge mis-id). The analogous pro¬ 
cesses p ± —» p ± e + e~ and p ± —> p ± p + p~ have negligible rates 
for the selected events. The tt and Z/y* —> 7 + 7~+jets events 
that undergo this process contribute to 270r| lad in the ee and 
ej± categories. As electrons pass through more material at high 
I 77 I, the charge mis-id rate increases as well, and so the elec¬ 
tron I 77 I < 1.37 requirement significantly reduces the impact of 
this background. The charge mis-id rate due to track curvature 
mismeasurement for electrons and muons is negligible. 

The charge mis-id probability is determined by a maximum- 
likelihood fit using Z —» ee events reconstructed as same-sign 
and as opposite-sign pairs, as a function of electron // and p-\. 
This probability function is then applied to a sample of events 
passing the 270Th a d selection except that the lepton pair is re¬ 
quired to be opposite sign. The charge mis-id probability from 
the relatively low momentum Z daughters is extrapolated to 
higher pj using scaling functions extracted from Monte Carlo 
simulations. The dominant uncertainty is due to the statistical 
precision of the charge mis-id probability determination, and is 
» 40% in the signal regions. 

7.4. Non-prompt light leptons 

A significant background arises from leptons not produced 
in decays of electroweak bosons (non-prompt leptons), which 
can promote (for example) a single-lepton tt event into a 2/7)T had 
category or a dilepton tt event to the 37 or 271 T| lad categor¬ 
ies. These backgrounds in the signal regions are expected to 
be dominated by tt or single top quark production with leptons 
produced in decays of heavy-flavour hadrons. Production of tt 
with an additional photon which converts in the detector mater¬ 
ial is a subdominant contribution. With the tight object selec¬ 
tion requirements applied in this analysis, almost all reconstruc¬ 
ted electron and muon objects correspond to real electrons and 
muons; the fraction arising from incorrect particle identification 
is negligible. Estimates of these backgrounds are obtained from 
data. Each channel has a slightly different procedure, motivated 
by the specific event topology and the statistical power avail¬ 
able in the control regions. The methods are discussed below, 
and the expected non-prompt lepton contributions to the vari¬ 
ous categories are shown in Table 3. In the following, a tight 
lepton is a lepton that passes the nominal selection, a sideband 
lepton is defined as a lepton candidate which satisfies different 
criteria than the tight lepton selection (identification selection, 
isolation, or p T ), and (sideband) control regions either require 
one or more sideband leptons to replace a tight lepton in the 


signal region selection, or have the same lepton selection as the 
signal region but different jet requirements. 

7.4.1. 270rhad categories 

The non-prompt lepton yields in the signal regions are 
estimated by extrapolating from sideband control regions in 
data which are enriched in tt non-prompt contributions. For 
electrons, sideband objects are selected by inverting the elec¬ 
tron identification and isolation requirements; for muons the 
sideband objects have low transverse momentum, 6 < pj < 
10 GeV, but otherwise are selected the same way as nominal 
muons. Transfer factors are used to extrapolate from events 
with one tight and one sideband lepton, but which otherwise 
pass the signal region selections, to the signal regions with two 
tight leptons. These transfer factors are determined from ad¬ 
ditional data control regions (tight + sideband and two tight 
leptons) with lower jet multiplicity (1 < « |et < 3 for electrons, 
2 < fljet < 3 for muons). In all regions the expected contribution 
from processes producing prompt leptons is subtracted before 
extracting transfer factors or using the yields for extrapolation. 
For channels with electrons, the charge mis-id background is 
also subtracted, and a dilepton mass veto is applied in the con¬ 
trol regions to suppress contributions from Z —» e + e" decays. 
A cross-check on the muon estimate, using an extrapolation in 
muon isolation instead of muon p T , agrees well with the nom¬ 
inal procedure and provides additional confidence in the estim¬ 
ate. 

The systematic uncertainties on this procedure are estim¬ 
ated by checking a) its ability to successfully predict the non¬ 
prompt background in tt simulation and b) the stability of the 
prediction using data when the selection of the control regions 
is altered. For the former, different parton shower and fc-hadron 
decay models were checked, as was the result of removing the 
/^-tagged jet requirement. In addition, for electrons, the effects 
of relaxing the pseudorapidity requirement to \r]\ < 2.5 and of 
raising the p-\ threshold were studied. These checks show sta¬ 
bility at the 25-30% level, limited by the statistical precision 
of the simulations. The stability in data is checked by alter¬ 
ing the pj required for the ^-tagged jet, applying a require¬ 
ment on missing transverse momentum 3 £™ ss , extracting the 
transfer factors only from events with three jets, or (for muons) 
using 10-15 GeV muons as the sideband objects. This check 
shows stability of the predictions to 14% for muons and 19% 
for electrons. Additional systematic uncertainties in the predic¬ 
tion arise from the statistical uncertainties on the yields in the 
control regions and the subtraction of prompt and charge mis- 
id contributions. The overall uncertainties on the non-prompt 
yield prediction in any given category range from 32% to 52%, 
and correlations between the categories due to uncertainties in 
the transfer factors are included in the fit (see Section 9). 


3 This is calculated using calorimeter energy deposits, calibrated according 
to associated reconstructed physics objects, and also including the transverse 
momenta of reconstructed muons. 
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Table 3: Expected and observed yields in each channel. Uncertainties shown are the sum in quadrature of systematic uncertainties 
and Monte Carlo simulation statistical uncertainties. “Non-prompt” includes the misidentified Th a d background to the 1 (2t\ vm \ 
category. Rare processes (fZ, ttWW, triboson production, tilt, tH ) are not shown as a separate column but are included in the total 
expected background estimate. 


Category 

q mis-id 

Non-prompt 

ttW 

ttZ 

Diboson 

Expected bkg. 

ttH (p = 1) 

Observed 

ee + >5 j 

1.1 ±0.5 

2.3 ± 1.2 

1.4 ±0.4 

0.98 ± 0.26 

0.47 ± 0.29 

6.5 ± 1.8 

0.73 ± 0.14 

10 

eg + > 5 j 

0.85 ± 0.35 

6.7 ± 2.4 

4.8 ± 1.2 

2.1 ±0.5 

0.38 ± 0.30 

15 ±3 

2.13 ± 0.41 

22 

+ ^ 5 j 

- 

2.9 ± 1.4 

3.8 ± 0.9 

0.95 ± 0.25 

0.69 ± 0.39 

8.6 ± 2.2 

1.41 ± 0.28 

11 

ee + 4 j 

1.8 ±0.7 

3.4 ± 1.7 

2.0 ± 0.4 

0.75 ± 0.20 

0.74 ± 0.42 

9.1 ±2.1 

0.44 ± 0.06 

9 

eg + 4 j 

1.4 ±0.6 

12 ±4 

6.2 ± 1.0 

1.5 ±0.3 

1.9 ± 1.0 

24 ± 5 

1.16 ±0.14 

26 

m + 4 j 

- 

6.3 ± 2.6 

4.7 ± 0.9 

0.80 ± 0.22 

0.53 ± 0.30 

12.7 ± 2.9 

0.74 ±0.10 

20 

3 1 

- 

3.2 ± 0.7 

2.3 ± 0.7 

3.9 ±0.8 

0.86 ± 0.55 

11.4 ±2.3 

2.34 ± 0.35 

18 

2^1 Thad 

- 

0 4 +o.6 
_0.4 

0.38 ±0.12 

0.37 ± 0.08 

0.12 ± 0.11 

1.4 ±0.6 

0.47 ± 0.08 

1 

lf2T ha d 

- 

15 ±5 

0.17 ±0.06 

0.37 ± 0.09 

0.41 ± 0.42 

16 ±5 

0.68 ±0.13 

10 

4 t Z-enr. 

- 

< 10~ 3 

< 3 x 10~ 3 

0.43 ± 0.12 

0.05 ± 0.02 

0.55 ±0.15 

0.17 ± 0.02 

1 

4t Z-dep. 

- 

< 10~ 4 

< 10~ 3 

0.002 ± 0.002 

< 2 x 10~ 5 

0.007 ± 0.005 

0.025 ± 0.003 

0 


7.4.2. 3 i category 

Sideband leptons are defined by reversing the isolation re¬ 
quirement for electrons and muons and, for electrons, requiring 
that the candidate fail the tight electron identification discrimin¬ 
ant requirement of the analysis but pass a looser selection. The 
non-prompt lepton contribution in the signal region is estim¬ 
ated by extrapolating from data regions with two tight and one 
sideband lepton, using transfer factors estimated from Monte 
Carlo simulation. These events typically contain two prompt 
opposite-sign leptons and one non-prompt lepton, which neces¬ 
sarily must be of the same sign as one of the prompt leptons. 
Therefore the non-prompt lepton estimation procedure is ap¬ 
plied only to the two same-sign leptons. The simulation-derived 
transfer factor is validated in a region of lower jet multiplicity 
(2 < Uj et < 3 and exactly one /^-tagged jet). Good agreement is 
observed in this validation region between the prediction ( 11.8 
± 2.3) and the observed yield (9.8 + 4.9 events after prompt 
background subtraction). Systematic uncertainties in the pro¬ 
cedure are derived by studying the agreement between data and 
simulation in the variables used for the extrapolation, which is 
« 20% for both electrons and muons. Additional uncertainties 
arise from the statistical uncertainties on the yields in the con¬ 
trol regions and in the tt simulation. 

7.4.3. 2(TThad category 

Reconstructing two same-sign light leptons from tt produc¬ 
tion or similar sources requires that one of the light leptons 
is non-prompt or has its charge misidentified. In the 2HThad 
category, the charge mis-id contribution is negligible and the 
primary concern is non-prompt light leptons. Around half of 
the Th a d candidates in these events come from W —* tv de¬ 
cays, while the remainder arise from misidentified light-quark 
or gluon jets. Regardless of whether the Thad candidate is a fake, 
there is also a non-prompt light lepton. Due to this fact, side¬ 
bands in the light-lepton selection criteria are used, analogously 
to the 2£0Th a d and 71 categories. Since the ratio of real and fake 
Thad candidates is similar in the signal and all control regions, 
fake Thad candidates are not accounted for separately; the small 


variations in the ratio in the control regions are found to have 
negligible impact on the total estimate in the signal region. In 
order to maintain similar origin composition of the non-prompt 
leptons, the E \ isolation requirement is inverted, the p-\ isola¬ 
tion requirement is relaxed, and for electrons the identification 
criteria are also relaxed to a looser working point. The low jet 
multiplicity region 2 < n ]ct < 3 is used to determine a transfer 
factor from sideband to tight lepton selections. The expected 
non-prompt lepton yield in the signal region is obtained by us¬ 
ing this transfer factor to extrapolate from a control region with 
the same jet selection as the signal region but with one tight 
and one sideband light lepton. The procedure is validated by 
checking that it correctly reproduces the signal region yield ex¬ 
pected in tt simulations. The assigned systematic uncertainty 
(27%) is dominated by the statistical precision of this test. The 
overall uncertainty on the non-prompt background prediction is 
dominated by the limited statistics of the high jet multiplicity 
control region. 

7.4.4. 4 1 category 

The non-prompt lepton contribution in this category is ex¬ 
pected to be negligible and is estimated to be < 10 3 events 
in the Z-enriched sample and < 10 4 events in the Z-depleted 
sample. In both cases this represents < 2% of the total back¬ 
ground expectation. These estimates are obtained using the 
transfer factors from the 3t channel and appropriate control 
regions with two loose leptons and relaxed jet multiplicity re¬ 
quirements. 

7.5. Thad misidentification in the l^Thad category 

The nominal estimate for the fake Thad yield is derived from 
tt simulation. To obtain a sufficiently large sample size, fast 
simulation using parameterized calorimeter showers is used. At 
all preselection stages the simulation is found to give an accept¬ 
able description of the tt background, both in kinematic distri¬ 
butions and total yield. This estimate is cross-checked with the 
data-driven method described below. 



Of the two rhad candidates, one is opposite in sign to the 
light lepton (OS) and the other has the same sign (SS). The SS 
candidate is almost always a fake Th a d- while the light lepton 
is prompt and the OS Th a d candidate is often real (~ 30%). A 
sideband Th a d is defined as a candidate passing a loose identific¬ 
ation BDT selection but not the nominal tight one. Assuming 
the Th a d candidate fake probabilities are not correlated between 
jets identified as OS and SS candidates, control regions can be 
used to predict yields in the signal region. There are three con¬ 
trol regions, depending on whether only the OS, only the SS, or 
both the OS and SS Th a d candidates are sideband objects. The 
two regions with sideband OS Th a d candidates are used to ob¬ 
tain the transfer factor for the SS Th a d candidate, which is then 
applied to the region with a tight OS and sideband SS candid¬ 
ate to obtain the prediction for the signal region where both are 
tight. The transfer factor is measured as a function of the p-y, 
77 , number of tracks, and /;-tag discriminant value of the SS Th a d 
candidate. The data-driven method is cross-checked in tt sim¬ 
ulation and found to successfully predict the yields in the sig¬ 
nal region. The main limitation of this method is the statistical 
power of the control regions. 

The simulation-driven method is taken as the primary es¬ 
timate, as the validation of the method at preselection stages is 
more precise than the data-driven method due to larger event 
yield for the former. The comparison of the simulation- and 
data-driven techniques gives a 36% uncertainty in the predic¬ 
tion in the signal region, which is taken as the systematic un¬ 
certainty on the estimate. 

8. Other systematic uncertainties 

Systematic uncertainties not already discussed are summar¬ 
ized below. 

The uncertainty on the integrated luminosity is 2.8%. This 
uncertainty is derived from a calibration of the luminosity 
scale derived from beam-separation scans performed in Novem¬ 
ber 2012 , following the same methodology as that detailed in 
Ref. [94], 

Lepton reconstruction and identification uncertainties are 
obtained from Z —> tt, Z —> tty, T —* tl, and //i^ —* tt 
events [80-82]. Uncertainties on the detector response are as¬ 
sessed similarly to other ATLAS analyses. The modelling of 
the efficiency of the tight isolation requirements in simulation 
is explicitly checked as a function of the number of jets in the 
event. These corrections are found to be very small, with un¬ 
certainties limited by data statistics. 

The largest jet-related systematic uncertainty arises from 
the jet energy scale, in particular contributions from the in-situ 
calibration in data, the different response to quark and gluon 
jets, and the pileup subtraction. The impact of the ^-tagging ef¬ 
ficiency uncertainty on the signal strength p = cr, f - Hobs /cr„- ftSM 
at the best-fit value of p is A p Because only one (of 

typically two) fo-jets present in signal or ttV events is required 
to be tagged, the uncertainty on the ^-tagging efficiency (while 
included) does not have as large an effect in this analysis as it 
does in other ttH searches such as those targeting the H —> bb 
decay. 


The uncertainties on the inclusive ttH production cross sec¬ 
tion are discussed in Section 3. Additionally, the effects of PDF 
uncertainty, QCD scale choice, and parton shower algorithm on 
the signal acceptance in each analysis category are considered. 
The resulting relative uncertainties on the acceptance are 0.3- 
1.4% for PDF, 0.1-2.7% for scale choice, and 1.5-13% for par- 
ton shower algorithm. 

For most backgrounds the uncertainties from Monte Carlo 
simulation sample sizes are negligible. For the diboson back¬ 
grounds, however, these can reach 50% of the total diboson 
yield uncertainties shown in Table 3. 

9. Results 

The observed yields, and a comparison with the expected 
backgrounds and ttH signal, are shown in Table 3. The distri¬ 
butions of the number of jets in the events passing signal re¬ 
gion selections are shown in Fig. 2. The best-fit value of the 
signal strength p = cr„- ffi 0 bs /cr,f^sM is determined using a max¬ 
imum likelihood fit to the data yields of the categories listed in 
Table 3, which are treated as independent Poisson terms in the 
likelihood. The fit is based on the profile-likelihood approach 
where the systematic uncertainties are treated as nuisance para¬ 
meters with prior uncertainties that can be further constrained 
by the fit [95], The p — 1 hypothesis assumes Standard Model 
Higgs boson production and decay with mu = 125 GeV; for 
all other values of p only the ttH production cross section is 
scaled (the Higgs boson branching fractions are fixed to their 
SM values). 

Systematic uncertainties are allowed to float in the fit as 
nuisance parameters and take on their best-fit values. The only 
constraints on nuisance parameter uncertainties found by the 
fit are for non-prompt lepton transfer factors and normaliza¬ 
tion region yields in the 2/'0r had categories and the fake T bad 
background yield in the 1 category. The former all have 

large statistical components and so the additional information 
from the signal regions is expected to constrain them. The lat¬ 
ter has a very large initial uncertainty which the fit is able to 
constrain as p is required to be the same in all categories. The 
largest difference between pre- and post-fit nuisance parameter 
values is in the 1^2rh a d fake estimate, which shifts by -l.Ocr 
due to the deficit of observed relative to expected events. The 
next largest effect is a +0.4cr shift in the 2({)t\ vm \ non-prompt p 
transfer factor. 

The results of the fit are shown in Fig. 3. The impact of 
the most important systematic uncertainties on the measured 
value of p in the combined fit is shown in Table 4. In each cat¬ 
egory, the uncertainties on p are mainly statistical, except for 
the combined 2 f 0 r had result where the statistical and system¬ 
atic uncertainties are of comparable size. In the At Z-depleted 
category, a (non-physical) signal strength p < -0.17 results in 
a negative expected total yield and a discontinuity in the pro¬ 
filed likelihood; the error bar is therefore truncated at this point. 
The results are compatible with the Standard Model expectation 
and with previous searches for ttH production in multilepton fi¬ 
nal states [18]. Combined over all categories, the value of p is 
found to be 2. In the presence of a signal of SM strength. 
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Figure 2: The spectrum of the number of jets expected and observed in each signal region. For display purposes the six 2fOT ha d 
categories (ee/e/V/#' and = 4/> 5 jets) are combined into one plot, as are the two 4f: categories (Z-enriched and Z-depleted). The 
hatched bands show the total uncertainty on the background prediction in each bin. The non-prompt and charge mis-id background 
spectra are taken from simulation of tt, single top, Z —> f + f“+jets, and other small backgrounds, with normalization as described 
in the text (in particular the = 4/> 5 jet regions of the 2{{)t\ vm \ pldOhave the ratio given by the data-driven prediction). The overlaid 
red line shows the ttH signal from the Standard Model. For visibility, the ttH signal is multiplied by a factor of 2.4 in the 
3£, and 1 C2Th a d plots. 





















































Table 5: Observed and expected 95% CL upper limits, derived using the CL S method, on the strength parameter p — cr ,tu,obs I c ttH ,sm 
for a Higgs boson of mass mu = 125 GeV. The last column shows the median expected limit in the presence of a ttH signal of 
Standard Model strength. 


Channel 

Observed Limit 

—2 cr 

—1 cr 

Expected Limit 
Median +lcr +2cr 

Median (p = 1) 

2 f0r h ad 

6.7 

2.1 

2.8 

3.9 

5.7 

8.4 

5.0 

3 £ 

6.8 

2.0 

2.7 

3.8 

5.7 

8.5 

5.1 

2 f’lr ha d 

7.5 

4.5 

6.1 

8.4 

13 

21 

10 

At 

18 

8.0 

11 

15 

23 

39 

17 

lf2r h ad 

13 

10 

13 

18 

26 

40 

19 

Combined 

4.7 

1.3 

1.8 

2.4 

3.6 

5.3 

3.7 


Table 4: Leading sources of systematic uncertainty and their 
impact on the measured value of p. 


Source 

Ap 

2 f0r ha d non-prompt muon transfer factor 

+0.38 

-0.35 

ttW acceptance 

+0.26 

-0.21 

ttH inclusive cross section 

+0.28 

-0.15 

Jet energy scale 

+0.24 

-0.18 

2/X)rhad non-prompt electron transfer factor 

+0.26 

-0.16 

ttH acceptance 

+0.22 

-0.15 

ttZ inclusive cross section 

+0.19 

-0.17 

ttW inclusive cross section 

+0.18 

-0.15 

Muon isolation efficiency 

+0.19 

-0.14 

Luminosity 

+0.18 

-0.14 


the combined fit is expected to return p = 1.0' [ ^ The p = 0 
hypothesis has an observed (expected) p-value of 0.037 (0.18), 
corresponding to 1.8cr (0.9cr). The p = 1 hypothesis (the SM) 
has an observed p-value of 0.18, corresponding to 0.9cr. The 
likelihood function can be used to obtain 95% confidence level 
(CL) upper limits on p using the C’L V method [95, 96], leading 
to the results in Table 5. The observed (expected) upper limit, 
combining all channels, is p < 4.7 (2.4). 

This analysis is a search for ttH production; as such, pro¬ 
duction of tHqb and tHW is considered as a background and set 
to Standard Model expectation. Including this contribution as a 
background induces a shift of Ap = -0.04 compared to setting 
it to zero. A full extraction of limits on the top quark Yukawa 
coupling including the relevant modifications of single top plus 
Higgs boson production is reported in Ref. [97]. 

The results are sensitive to the assumed cross sections for 
ttW and ttZ production, and use theoretical predictions for these 
values as experimental measurements do not yet have sufficient 
precision. The best-fit p value as a function of these cross sec¬ 
tions is 
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Figure 3: Best-fit values of the signal strength parameter 
P = cr„- Hi0bs /cr„- HSM . For the At Z-depleted category, p < -0.17 
results in a negative expected total yield and so the lower un¬ 
certainty is truncated at this point. 
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10. Conclusions 

A search for ttH production in multilepton final states has 
been performed using 20.3 fb 1 of proton-proton collision data 
at = 8 TeV recorded by the ATLAS experiment at the LHC. 
The best-fit value of the ratio fi of the observed production rate 
to that predicted by the Standard Model is This res¬ 

ult is consistent with the Standard Model expectation. A 95% 
confidence level limit of f t < 4.7 is set. The expected limit in 
the absence of ttH signal is fi < 2.4. The observed (expected) 
/;-value of the no-signal hypothesis corresponds to 1.8cr (0.9cr). 
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